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Q^. Abstract 

■ We study the regularity of the extremal solution of the semilinear biharmonic equation A'^u = ^y^^i 

lO i which models a simple Micro-Electromechanical System (MEMS) device on a ball B C R^, under 

Dirichlet boundary conditions u = dvU = on dB. We complete here the results of F.H. Lin and Y.S. 
Yang [M] regarding the identification of a "pull-in voltage" A* > such that a stable classical solution 
', Ma with < u\ < 1 exists for AG (0, A*), while there is none of any kind when A > A*. Our main result 

• asserts that the extremal solution u\* is regular (sup^ u\* < 1) provided N < 8 while u\* is singular 

' (sup^ u\* — 1) for A*' > 9, in which case x]"^^"^ on the unit ball, where 

'S: Co:= (^)* andA:=|(7V-|)(iV-|). 



1 Introduction 

The following model has been proposed for the description of the steady-state of a simple Electrostatic 
^ ! MEMS device: 

, ^ < u < 1 inn '^^f 

'nI" ' Ku ~ adyU — on dQ, 

o ■ 

' where 7,x > 0, / G C(0, [0,1]) are fixed, O is a bounded domain in and A > is a varying 

parameter (see for example Bernstein and Pelesko [13). The function u{x) denotes the height above a point 
X € n C 'R^ of a dielectric membrane clamped on dfl, once it deflects forwards a ground plate fixed at 
, height z — I, whenever a positive voltage - proportional to A - is applied. 

■ In studying this problem, one typically makes various simplifying assumptions on the parameters a,/3, 7,x, 

and the first approximation of ^ that has been studied extensively so far is the equation 

r 1J.1 a u 
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where we have set a = /? = x = and 7=1 (see for example [51 [HI IS] and the monograph [71) . This simple 
model, which lends itself to the vast literature on second order semilinear eigenvalue problems, is already 
a rich source of interesting mathematical problems. The case when the "permittivity profile" / is constant 
(/ = 1) on a general domain was studied in [TBI, following the pioneering work of Joseph and Lundgren 
[T3] who had considered the radially symmetric case. The case for a non constant permittivity profile / was 
advocated by Pelesko [18,^ , taken up by [TT] , and studied in depth in [HI HI ISl • The starting point of the 
analysis is the existence of a pull-in voltage A* (17, /), defined as 



It is then shown that for every < A < A*, there exists a smooth minimal (smallest) solution of (<5')a,/, while 
for A > A* there is no solution even in a weak sense. Moreover, the branch A i— > ux{x) is increasing for each 
a; G il, and therefore the function u*{x) := Hm^/A* "^(a;) can be considered as a generalized solution that 
corresponds to the pull-in voltage A*. Now the issue of the regularity of this extremal solution - which, by 
elliptic regularity theory, is equivalent to whether sup^ m* < 1 - is an important question for many reasons, 
not the least of which being the fact that it decides whether the set of solutions stops there, or whether a new 
branch of solutions emanates from a bifurcation state (u*,A*). This issue turned out to depend closely on 
the dimension and on the permittivity profile /. Indeed, it was shown in [9] that u* is regular in dimensions 
1 < < 7, while it is not necessarily the case for N > 8. In other words, the dimension = 7 is critical 
for equation {S)\ (when / = 1, we simplify the notation (5')a,i into {S)\). On the other hand, it is shown 
in [5] that the regularity of u* can be restored in any dimension, provided we allow for a power law profile 
a; I'' with r] large enough. 

The case where /3 = 7 = x = (and a = 1) in the above model, that is when we are dealing with the 
following fourth order analog of {S)\ 



was also considered by [H [T31 but with limited success. One of the reasons is the lack of a "maximum 
principle" which plays such a crucial role in developing the theory for the Laplacian. Indeed, it is a well 
known fact that such a principle does not normally hold for general domains Q (at least for the clamped 
boundary conditions u — d^u = on dO) unless one restricts attention to the unit ball Q, = B m , 
where one can exploit a positivity preserving property of due to T. Boggio [3 . This is precisely what 
was done in the references mentioned above, where a theory of the minimal branch associated with {P)\ is 
developed along the same lines as for {S)\. The second obstacle is the well-known difhculty of extracting 
energy estimates for solutions of fourth order problems from their stability properties. This means that the 
methods used to analyze the regularity of the extremal solution for {S)\ could not carry to the corresponding 
problem for {P)\. 

This is the question we address in this paper as we eventually show the following result. 

Theorem 1.1. The unique extremal solution u* for {P)\* in B is regular in dimension 1 < N < 8, while it 
is singular (i.e, sup^ u* = 1) for N > 9. 

In other words, the critical dimension for {P)\ in B is N = 8, as opposed to being equal to 7 in {S)\. We 
add that our methods are heavily inspired by the recent paper of Davila et al. [5] where it is shown that 
= 12 is the critical dimension for the fourth order nonlinear eigenvalue problem 



while the critical dimension for its second order counterpart (i.e., the Gelfand problem) is = 9. There is, 
however, one major difference between our approach and the one used by Davila et al. [5j. It is related to the 



A*(r2, /) := sup A > : there exists a classical solution of {S)xj 





in B 
on dB, 
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most delicate dimensions - just above the critical one ~ where they use a computer assisted proof to establish 
the singularity of the extremal solution, while our method is more analytical and relies on improved and non 
standard Hardy-Rellich inequalities recently established by Ghoussoub-Moradifam [TD] (See Appendix). 

Throughout this paper, we will always consider problem {P)\ on the unit ball B. We start by recalling some 
of the results from 4] concerning {P)\, that will be needed in the sequel. We define 

A* :— sup |a > : there exists a classical solution of (P)a|, 

and note that we are not restricting our attention to radial solutions. We will deal also with weak solutions: 
Definition 1.1. We say that u is a weak solution of {P)\ if < u < 1 a.e. in B, (^i^^yi G and 

uA'^ - A / V0 e C^B) n H^{B). 

Jb U ~ 

We say that u is a weak super- solution (resp. weak sub-solution) of {P)\ if the equality is replaced with the 
inequality > (resp. <) for all (j) G C"^(B) n Hq{B) with (j)>0. 

We also introduce notions of regularity and stability. 

Definition 1.2. Say that a weak solution u of {P)\ is regular (resp. singular) if \\u\\oo < 1 (resp. =) and 
stable (resp. semi-stable) if 

fi.iu) = inf I J^iA^f - 2XJ^^^^-^ : e H^.iB), ||0|U. = l| 

is positive (resp. non-negative). 

The following extension of Boggio's principle will be frequently used in the sequel (see [51 Lemma 16] and 
[3 Lemma 2.4]): 

Lemma 1.1 (Boggio's Principle). Let u G L^{B). Then u > a.e. in B, provided one of the following 
conditions hold: 

1. ue C^(B), A^u >0 onB, andu^ §^ = on OB. 

2. Jg uA'^(j)dx >0 for alio <(j)e C^iB) n H^{B). 

3. ue H^{B), u = and^^<Q on dB, and AuA0 > for all < (j) e H^iB). 
Moreover, either u = or u > a.e. in B. 

The following theorem summarizes the main results in W that will be needed in the sequel: 
Theorem 1.2. The following assertions hold: 

1. For each < A < A* there exists a classical minimal solution u\ of {P)\. Moreover u\ is radial and 
radially decreasing. 

2. For A > A*, there are no weak solutions of {P)\. 

3. For each x e B the map A ^ ux{x) is strictly increasing on (0, A*). 
4- The pull-in voltage A* satisfies the following bounds: 

J ^ i -f Af \f AT <r \* <r i 



s ,-(Af- -)(7V- -) ^ < A* < 

\ 27 9^ 3^^ 3V ~ ~ 27 



where v\ denotes the first eigenvalue of A^ in Hq (B) . 
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5. For each < A < A*, u\ is a stable solution (i.e., ^i{u\) > Q). 

Using the stability of ux, it can be shown that u\ is uniformly bounded in Hq{B) and that is uniformly 
bounded in L^{B). Since now A i— > u\{x) is increasing, the function u*{x) := limA/^A* ux{x) is well defined 
(in the pointwise sense), u* G Hq{B), j-t— t G L^[B) and u* is a weak solution of {P)\*. Moreover u* is the 
unique weak solution of {P)\*. 

The second result we list from [4 is critical in identifying the extremal solution. 

Theorem 1.3. // u G Hl(B) is a singular weak solution of {P)x, then u is semi- stable if and only if 
(w,A) = (w*, A*). 



2 The effect of boundary conditions on the pull-in voltage 

As in [5j, we are led to examine problem {P)\ with non- homogeneous boundary conditions such as 

A^u = 7t^4t in B 

a < u < 1 in B {P)\,a,f} 

u ^ a , duU ~ (3 on dB, 



where a, (3 are given 

Notice first that some restrictions on a and /3 are necessary. Indeed, letting $(a;) := (a — ^) + ^\x\'^ denote 



the unique solution of 

A2$ = in B 



$ = a , 9^$ = /3 on dB, "-^^ 
we infer immediately from Lemma 1 1.11 that the function u — $ is positive in B, which yields to 

sup $ < sup u < 1 . 

B B 

To insure that <i> is a classical sub-solution of {P)\.a,f3, we impose a ^ 1 and (3 < 0, and condition sup$ < 1 

B 

rewrites as a — ^ < 1. We will then say that the pair (a, /3) is admissible if /? < 0, and a — ^ < 1. 

This section will be devoted to obtaining results for {P)\,a,f3 when (a,/3) is an admissible pair, which are 
analogous to those for {P)\. To cut down on notation, we shall sometimes drop a and (3 from our expressions 
whenever such an emphasis is not needed. For example in this section u\ and u* will denote the minimal 
and extremal solution of {P)\,a,i3- 

We now introduce a notion of weak solution for (-P)a,q./3- 

Definition 2.1. We say that u is a weak solution of {P)x^a.f3 if a < u < 1 a.e. in B, ^iJ^^y^ ^ ^^{^) ^'^^ 

f {u- $)A20 = A / V0 G C\B) n H^{B), 

Jb JB[^-uy 

where $ is given in 0). We say that u is a weak super- solution (resp. weak sub-solution) of {P)\^a.[3 if the 
equality is replaced with the inequality > (resp. <) for 4> >0. 

We now define as before 

A* sup{A > : {P)\.a,p has a classical solution} 

and 

:~ sup{A > : {P)\.a,fj has a weak solution}. 

Observe that by the Implicit Function Theorem, one can always solve {P)\^a,f) for small A's. Therefore, A* 
(and also A*) is well defined. 

Let now [/ be a weak super-solution of {P)\,a,f}- Recall the following standard existence result. 
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Theorem 2.1 ([2]). For every < / G L^{B), there exists a unique < u £ L^{B) which satisfies 

Jb Jb 
for all(l)eC*(B)nH^){B). 

We can now introduce the following "weak iterative scheme": Start with uq ~ U and (inductively) let u„, 
71 > 1, be the solution of 

/ K - = A / V (A e C\B) n Hl{B) 

Jb Jb (1 - Un-iY 

given by Theorem l2.1l Since is a sub-solution of (P)a,q,/3, one can easily show inductively by using Lemma 
11.11 that a < Un+i < Un < U for every n >0. Since 

(l-u„)-2 < (i-u)-^ eL\B), 

we get by Lebesgue Theorem, that the function u lim u„ is a weak solution of {P)\.a,f3 such that 
a < u <U . In other words, the following result holds. 

Proposition 2.1. Assume the existence of a weak super-solution U of {P)\,a,i3- Then there exists a weak 
solution u of {P)\.a,i3 so that a < u < U a.e. in B. 

In particular, we can find a weak solution of {P)\^a.i3 for every A G (0, A,). Now we show that this is still 
true for regular weak solutions. 

Proposition 2.2. Let {a, (3) he an admissible pair and let u be a weak solution of {P)\^a,f3- Then for every 
< < A, there is a regular solution for {P)fj,,a,i3- 

Proof. Let e G (0, 1) be given and let u — {1 — e)u + e$, where $ is given in ([2]). We have that 
supu <(!—£) + esup^ < 1 , inf M > (1 — e)a + einf $ = a, 

B B B B 

and for every < G C"'(i?) n Hl{B) there holds: 

[{u-^)A^(l) = (1 - e) /" (u - $)A20 = (1 - e)A /" 
Jb Jb Jb 
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{1~U) 

(l-^)'A/ > (1 - e)3A 



(1 -u + e(4'- 1))2 - ' ' Jsil-uf' 

Note that < (1 — e)(l — u) = 1 — m + £(<!> — 1) < 1 — m. So m is a weak super-solution of {P){i-e)^\.a,i3 
satisfying supw < 1. From Proposition 12.11 we get the existence of a weak solution w of {P)ii-ey^\ a p so 

B 

that a < w < u. In particular, supw < 1 and w is a regular weak solution. Since e G (0, 1) is arbitrarily 

B 

chosen, the proof is complete. □ 



Proposition 12.21 implies in particular the existence of a regular weak solution Ux for every A G (0, A*). 
Introduce now a "classical" iterative scheme: uq = and (inductively) u„ — Vn'V'^, n>l, where u„ G Hq{B) 
is the (radial) solution of 

A2i,„ = A2(w„-$) = inB. (3) 

(1 - Un^iY 

Since w„ G Hq{B), Un is also a weak solution of ([3]), and by Lemma Fl. II we know that a <Un < Wn+i < U\ 
for every n > 0. Since supu„ < supC/A < 1 for n > 0, we get that (1 — u„_i)^^ G L'^{B) and the 
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existence of w„ is guaranteed. Since Vn is easily seen to be uniformly bounded in H'^{B), we have that 
u\ :— lim Un does hold pointwise and weakly in H^{B). By Lebesgue Theorem, we have that ux is a 

n — >+OD 

radial weak solution of {P)x,a.f:i so that supwA < sup [/a < 1. By elliptic regularity theory fT] ux £ C°°{B) 

B b 

and Ma — $ = d,y{ux — $) = on dB. So we can integrate by parts to get 

/ A^ux(f>= [ A\ux - <i>)(f> ^ [ {ux-^)A^c^ = X [ . ^ „ 
Jb Jb Jb Jb (1 - uxY 

for every cj) G C^{B) n Hq{B). Hence, ux is a radial classical solution of {P)x.a,p showing that A* = A*. 
Moreover, since $ and vx := ux — ^ are radially decreasing in view of [201, we get that ma is radially 
decreasing too. Since the argument above shows that ux <U for any other classical solution U of {P)^,a,i3 
with > A, we have that ux is exactly the minimal solution and ux is strictly increasing as A t A*. In 
particular, we can define u* in the usual way: u*{x) — lim ux{x). 

Finally, we show the finitcncss of the pull-in voltage. 

Lemma 2.1. If {a, (3) is an admissible pair, then X*{a,f3) < +oo. 

Proof. Let M be a classical solution of {P)x,a,i3 and let {ifj, vi) denote the first eigenpair of in Hq{B) with 
ip > 0. Now, let C be such that 

/ {f3A'tP - ad^Aij:) = C f ip. 
Job Jb 

Multiplying {P)x.a,p by ip and then integrating by parts one arrives at 

/.((T^-"--^)* = °- 

Since ip > there must exist a point x & B where (-x_J(e))2 ~ i'iu{x) — C < 0. Since a < u{x) < 1, one can 
conclude that A < suPq,<„<]^(i^iu + C)(1 — u)^, which shows that A* < +cxd. □ 

The following summarizes what we have shown so far. 

Theorem 2.2. // {a, f3) is an admissible pair, then A* := A*(a,/3) € (0, +oo) and the following hold: 

1. For each < A < A* there exists a classical, minimal solution ux of {P)x,a,p- Moreover ux is radial 
and radially decreasing. 

2. For each x d B the map A t— > ux{x) is strictly increasing on (0, A*). 

3. For A > A* there are no weak solutions of {P)x.a,i3- 



2.1 Stability of the minimal branch of solutions 

This section is devoted to the proof of the following stability result for minimal solutions. We shall need 
yet another notion of H^{B)~weak solutions, which is an intermediate class between classical and weak 
solutions. 

Definition 2.2. We say that u is a H^{B) — weak solution of (P)a,q,/3 if u — <^ £ H^{B), a < u < 1 a.e. in 

B, j^£L\B) and if 

[ AuAcj, = X f ^ V0 G C\B) n H^{B), 

Jb Jb U ~ 

where $ is given in 0). We say that u is a H^{B)~weak super- solution (resp. H^(B) — weak sub-solution) 
iP)x.a,p if for (/) > the equality is replaced with > (resp. <) and u > a (resp. <), d^u < (3 (resp. >) 
on dB. 
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Theorem 2.3. Suppose (a,/3) is an admissible pair. 

1. The minimal solution u\ is then stable and is the unique semi-stable H^{B)~weak solution of {P)x^ci,f}- 

2. The function u* := lim u\ is a well-defined semi-stable H^{B) — weak solution of (P)a*,q,/3- 

3. When u* is classical solution, then /ii(it*) ~ and u* is the unique H^{B) — weak solution of {P)\* ^a,0- 
4- If V is a sinqular, semi-stable H^{B) — weak solution of {P)\^a.p- then v — u* and A = A* 

The crucial tool is a comparison result which is valid exactly in this class of solutions. 

Lemma 2.2. Let (a,/3) be an admissible pair and u be a semi-stable H^{B) — weak solution of (P)a,q./3- 
Assume U is a H^(B) — weak super- solution of {P)\.a,[3 so that J7 — $ G Hq{B). Then 

1. u < U a.e. in B; 

2. If u is a classical solution and fJ.i{u) = then U = u. 

Proof, (i) Define w := u — U. Then by the Moreau decomposition [T7] for the biharmonic operator, there exist 
wi,W2 e Hq{B), with w = wi -\- wi > Q a.e., A^W2 < in the iif^(_B)— weak sense and Jg AwiAw2 = 0. 
By Lemma ri.li we have that i/;2 < a.e. in B. 
Given now < (/> £ C^{B), we have that 

/ AwAc^<X f {f{u)-f{U))c^, 

J B J B 

where f{u) := (1 — Since u is semi-stable, one has 

A / f\u)wl < f {Awif = / AwAwi < A / {f{u)-f{U))wi. 



Since wi > w one also has 



which once re-arranged gives 



jj'{u)ww, < I {f{u)~f{U))Wi, 



fwi > 0, 



where f{u) — f{u) — f{U) — f'{u){u — U). The strict convexity of / gives / < and / < whenever u . 
Since wi>Q a.e. in B one sees that w <Q a.e. in B. The inequality u < U a.e. in B is then established. 



ii) Since m is a classical solution, it is easy to see that the infimum in is attained at some 0. The 

2A 

(l-tl)3 



function (p is then the first eigenfunction of — ^ ,3 in Hq{B). Now we show that is of fixed sign. 



Using the above decomposition, one has = 01 -f 02 where (pi G IIq{B) for i = 1, 2, 0i > 0, A(j)iA(j)2 = 
and A'^(j)2 < in the iJQ(i3)— weak sense. If (j) changes sign, then 01 ^ and 02 < in B (recall that either 
02 < or 02 = a.e. in B). We can write now: 

^ ^ ^ /^(A(0i - 02))^ - Af H(0i - 02)^ ^ ^(A0)^ - Xf{u)<p^ 
Jb(0i-02)^ Jb0^ 

in view of 0i02 < — 0i02 in a set of positive measure, leading to a contradiction. 

So we can assume > 0, and by the Boggi's principle we have > in _B. For < t < 1 define 

g{t)^ I A[tU +[l-t)u]A<f)-\ ( f{tU +[l~t)u)(j), 
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where 4> is the above first eigenfunction. Since / is convex one sees that 

g{t) > A / [tf{U) + (1 - t)fiu) - fitu + (1 - t)u)] > 
Jb 

for every t > 0. Since .g(0) = and 

g'(0)= / AiU -u)A^-Xf'{u){U -u)(f> = 0, 



we get that 

<?"(0) = -A / f"{u){U-uf^>0. 



Since f"{u)4> > in i?, we finaUy get that U = u a.e. in B. □ 
Based again on Lemma [1. IT S), we can show a more general version of the above Lemma [2T2l 



Lemma 2.3. Let (a,/3) be an admissible pair and j3' < 0. Let u be a semi-stable H^{B)~weak sub-solution 
of {P)x^a.i3 with u = a, d^u = (3' > (3 on dB. Assume that U is a H^(B) — weak super- solution of {P)\,a,p 
with U = a, d^U = P on dB. Then U > u a.e. in B. 

Proof. Let u G Hq{B) denote a weak solution to A^w = A^(m — U) in B. Since u — u + U = and 
9i/({t — M + [/) < on dB, by Lemma 11.11 one has that u > u — U a.e. in B. Again by the Moreau 
decomposition [17] . we may write u as u = w -\- v, where w,v £ Hq{B), w > a.e. in B, A^w < in a 
ij2(B)-weak sense and /g AwAv = 0. Then for < G C"'(B) n H^{B) one has 



AuA(l)^ / A{u~U)A(t) < X / {f{u)^f{U))(l). 
Jb Jb 

f AuAw <X f {f{u) - f{U))w. 
Jb Jb 

X [ f{u)w^ < [ (Aw)2 = / AuAw, 
Jb Jb Jb 



In particular, we have that 
Since by semi-stability of u 

we get that 

/ f'{u)w' < f {f{u)-f{U))w. 

J B J B 

By Lemma ll.ll we have v <0 and then w > u > u — U a.e. in B. So we see that 

0< / if{u)~fiU)~f'iu)iu-U))w. 



The strict convexity of / implies as in Lemma 12^ that U > u a.e. in B. □ 
We shall need the following a-priori estimates along the minimal branch ux. 
Lemma 2.4. Let {a, (3) be an admissible pair. Then one has 

{ux ^ [ ux-^ 



{1-uxf - Jb {l~uxr 
where $ is given in In particular, there is a constant C > so that for every X G (0, A*), we have 
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Proof. Testing {P)\,a,f3 on ua - $ G C^{B) n H^{B), we see that 

^/ Tr^= f A^aAK - $) = / (AK-$))^>2A / ^^^^^ 
in view of = 0. In particular, for S > small we have that 



1 < 1 / 1 /■ 1 



J {\u^-<i.\>5} (1 - 

by means of Young's inequality. Since for 5 small, 

1 



<C' 



'{I«A-*|<5} (1 ^ '"^)^ 

for some C > 0, we can deduce that for every A G (0, A*), 
for some C > 0. By Young's and Holder's inequalities, we now have 



and estimate ^ is therefore established. □ 
We are now ready to establish Theorem 12.31 

Proof (of Theorem 12. 3p : (1) Since ||ua||oo < 1, the infimum defining ^i{ux) is achieved at a first 
eigenfunction for every A G (0, A*). Since A ^ ux{x) is increasing for every a; G it is easily seen that 
A ^ fJ-i{ux) is an increasing, continuous function on (0, A*). Define 

A,» := sup{0 < A < A* : Hi{ux) > 0}. 

We have that A*» — A*. Indeed, otherwise we would have that fii{ux,,) = 0, and for every fj> G (A^f:^, A*) 
would be a classical super-solution of {P)x,,,a.p- A contradiction arises since Lemma 12.21 implies — ux,,. 
Finally, Lemma [2?2l guarantees uniqueness in the class of semi-stable 7J^(_B)— weak solutions. 

(2) By estimate ^ it follows that ux — > u* in a pointwise sense and weakly in H^{B), and j-r^ G L^{B). 
In particular, u* is a iJ^(_B)— weak solution of (P)a'.q,/3 which is also semi-stable as limiting function of the 
semi-stable solutions {ux}- 

(3) Whenever ||w*||oo < 1, the function u* is a classical solution, and by the Implicit Function Theorem we 
have that ^i{u*) = to prevent the continuation of the minimal branch beyond A*. By Lemma [2.21 u* is 
then the unique iJ'^(_B)— weak solution of {P)x',a,(3- An alternative approach -which we do not pursue here- 
based on the very definition of the extremal solution u* is available in ^ when a = P = (see also (TS]) to 
show that u* is the unique weak solution of (P)a», regardless of whether u* is regular or not. 

(4) If A < A*, by uniqueness v = ux- So u is not singular and a contradiction arises. 

By Theorcm l2.2f 3') we have that A = A*. Since is a semi-stable iJ^(_B)— weak solution of (-P)a*,q,/3 and u* 
is a iJ^(_B)— weak super-solution of {P)x-',a,i3, we can apply Lemma [2.21 to get v < u* a.e. in B. Since u* is 
a semi-stable solution too, we can reverse the roles of v and u* in Lemma [2.21 to see that v > u* a.e. in B. 
So equality v — u* holds and the proof is done. 
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3 Regularity of the extremal solution for 1 < N < 8 

We now return to the issue of the regularity of the extremal solution in problem (-P)a- Unless stated 
otherwise, ux and u* refer to the minimal and extremal solutions of {P)\. We shall show that the extremal 
solution u* is regular provided 1 < < 8. We first begin by showing that it is indeed the case in small 
dimensions: 

Theorem 3.1. u* is regular in dimensions 1 < N < 4. 

Proof. As already observed, estimate ^ implies that f{u*) = (1 — u*)^^ G Li[B). Since u* is radial and 
radially decreasing, we need to show that u*(0) < 1 to get the regularity of u* . The intcgrability of f{u*) 
along with elliptic regularity theory shows that u* G iy*'2 (_B). By the Sobolev imbedding Theorem we get 
that u* is a Lipschitz function in B. 
Now suppose u*(0) = 1 and 1 < < 3. Since 

1 C 

> — in B 



1 ~ u \x\ 

for some C > 0, one sees that 



A contradiction arises and hence u* is regular for 1 < < 3. 

For iV = 4 we need to be more careful and observe that u* G (i?) by the Sobolev Imbedding Theorem. 
If u*(0) = 1, then Vm*(0) = and 

1 C . ^ 
>—^ mB 

l~U- \x\a 

for some C > 0. We now obtain a contradiction exactly as above. □ 
We now tackle the regularity of u* for 5 < < 8. We start with the following crucial result: 
Theorem 3.2. Let N > 5 and (u* , A*) be the extremal pair of (P)a- When u* is singular, then 

l-u*(a;) < Co|a;|^ inB, 

where Co := (^)' and X = \n := |(Af - |)(A^ - |). 
Proof. First note that Theorem 1 1.2r 4') gives the lower bound: 

\*>\=liN^^)iN^^). (5) 

For (5 > 0, we define us{x) := 1 — C5|a;|3 with Cs (x^"*^)^ ^ ^i^^e A^ > 5, we have that 
«5 6 iJf„^(R^), G i?oc(I^^) and ug is a H^-weak solution of 

A'us = ^^ inR^. 
(1 - usY 

We claim that us < u* in B, which will finish the proof by just letting (5^0. 

Assume by contradiction that the set F :— {r G (0,1) : us{r) > u*{r)} is non-empty, and let ri — sup F. 
Since 

usil)^l-Cs <0 = u*{l), 



1 
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we have that < ri < 1 and one infers that 

a := u*(ri) = us{ri) , (3 (u*)'(ri) > u^(ri). 

_ 4 

Setting U5,ri (f") = ^ (^^(rir) — 1)+1, we easily see that U5,ri is a weak super-solution of 

where ^ ^ 

a' := rj"^(a - 1) + 1 , P'-.^r^^p. 

_ 4 

Similarly, let us define u*^{r) = r^ ^ {u*{rir) — 1) + 1. The dilation map 

_ 4 

w — > Wri (r) — {w{rir) — 1) + 1 (6) 
is a correspondence between solutions of {P)\ on B and of (P) _4 on B -i which preserves the 

A,l— ^ ,0 1 

integrability. In particular, (m* , A*) is the extremal pair of (P) _4 on B^-i (defined in the 

A,l — ^ ,0 1 

obvious way). Moreover, u*_^ is a singular semi-stable H^{B)— weak solution of (P)A*,a',/3'- 

Since u* is radially decreasing, we have that 13' < 0. Define the function w as w{x) (a' — ^) + ^|a;|^+7(x), 
where 7 is a solution of A^7 = A* in i? with 7 = 9iy7 = on dB. Then u; is a classical solution of 

A'^w = A* in B 



J7 = a' , di,w ~ on (9i3. 

Since -^^j > A*, by Lemma [TTT] we have u*^ > w a.e. in _B. Since w{0) = a' — ^ + 7(0) and 7(0) > 0, 

the bound u*_^ < 1 a.e. in B yields to a' - < 1. Namely, (a',/3') is an admissible pair and by Theorem 
I2.3r 4) we get that {u*^,X*) coincides with the extremal pair of (P)A,a',/3' in B. 

Since (a', (3') is an admissible pair and ug.n is a weak super-solution of (P)a*+5Aiv q' /3'' Propo- 

sition [^TT] we get the existence of a weak solution of (P)a»+5Ajv a' p'- Since A* -t- SXn > A*, we contradict the 
fact that A* is the extremal parameter of (P)a,q',/3'- D 

Thanks to this lower estimate on m*, we get the following result. 

Theorem 3.3. //5 < iV < 8, then the extremal solution u* of {P)\ is regular. 

Proof. Assume that u* is singular. For e > set ^p{x) := and note that 

(AV')2 = (ffAr + 0(e))|xr^+2e 

where 

iV2(Ar-4)2 
Hn ■ 



16 

Given 77 £ C^{B), and since iV > 5, we can use the test function rjtp G Hq{B) into the stability inequality 
to obtain 



where 0(1) is a bounded function as e \ 0. By Theorem 13.21 we find that 

2~Xn f ^< f (A^r + 0(1), 
Jb fI Jb 

and then 

2AAr / \x\-^+'-' <{HM + 0{e)) f \x\-^+^' + 0(1). 



Computing the integrals one arrives at 

2Ajv < HN + 0{e). 

As e — > finally we obtain 2AAr < H^. Graphing this relation one sees that TV > 9. □ 
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We can now slightly improve the lower bound 
Corollary 3.1. In any dimension N > 1, we have 

\*>\N^liN-^)iN-^). (7) 

Proof. The function u :— 1 ~ \x\^ is a H^{B)— weak solution of {P)\j^ o -|- by contradiction A* ~ Xn, 
then u is a H^{B)—weak super-solution of {P)\ for every A G (0, A*). By Lemma [2^ we get that u\ < u for 
ah A < A*, and then u* < u a.e. in B. 

If 1 < TV < 8, u* is then regular by Theorems [XT] and By Theorem (23); 3) there holds fii{u*) = 0. 
Lemma 12.21 then yields that u* — u, which is a contradiction since then u* will not satisfy the boundary 
conditions. 

If now iV > 9 and A = A*, then Co = 1 in Theorem 13. 2( and we then have u* > u. It means again that 
u* = u, a contradiction that completes the proof. □ 



4 The extremal solution is singular for > 9 

We prove in this section that the extremal solution is singular for > 9. For that we have to distinguish 
between three different ranges for the dimension. For each range, we will need a suitable Hardy-Rellich type 
inequality that will be established in the appendix, by using the recent results of Ghoussoub-Moradifam [10] . 
As in the previous section (u*. A*) denotes the extremal pair of (P)a- 

• Case N > 17: To establish the singularity of u* for these dimensions we shall need the following well 
known improved Hardy-Rellich inequality, which is valid for N > 5. There exists C > 0, such that for all 

<t>eHi{B) 

(A0)2 dx > / ^dx + C [ cf'dx. (8) 



16 



B 



Case 10 < < 16: For this case, we shall need the following inequality valid for all (f) G Hq{B) 
^ ^ " 16 Jb i\x\^ -\x\^+'){\x\^ -\x\f) 



{N-l){N -A)^ 



2 r ja 



• Case = 9: This case is the trickiest and will require the following inequality for all G Hq{B) 

(A0f > ^ Q{\x\) {p{\x\) + ^) ^\ (10) 
where 

P(r) = and Q{r) ~ 



with if and ip being two appropriately chosen polynomials, namely 

^(r) := r-f +1 +r - 1.9 and i;{r) := r-^+^ + 20r-^-'^^ + lOr-'^ + lOr + 7r^ - A8. 

Recall that for a radial function tp, we set ANip{r) — ip"{r) + ^^^-^ip'{r). 

We shall first show the following upper bound on u*. 
Lemma 4.1. IfN>9, then u* < 1 - \x\^ in B. 
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Proof Recall from Corollary [O that A := |(iV - |)(iV - |) < A*. Wc now claim that ux < u for aU 
A 6 (A, A*). Indeed, fix such a A and assume by contradiction that 

Ri := inf{0 < R < 1 : ux < u in {R, 1)} > 0. 

From the boundary conditions, one has that ux{r) < u{r) as r ^ 1^. Hence, < i?i < 1, a := ux{Ri) — 
u{Ri) and f3 := u'x{Ri) < u'{Ri). Introduce, as in the proof of Theorem 13.21 the functions ux.Ri and urj^. 
We have that ux,Ri is a classical super-solution of (^')Ajv,ct',/3'' where 

a' := (a - 1) + 1 , /?' := R'^^p. 

Note that ur^ is a H^{B)—weak sub-solution of {P)x,^ a' p' which is also semi-stable in view of the Hardy- 
Rellich inequality ([5]) and the fact that 

2A <7T .^^'(^-4)^ 



16 

By Lemma [^31 we deduce that ux.Ri ^ ^Ri in B. Note that, arguing as in the proof of Theorem l3.21 (a', /3') 
is an admissible pair. We have therefore shown that ua > u in Br-^^ and a contradiction arises in view of the 
fact that lim u{x) — 1 and ||ua||oo < 1- It follows that ux < u in B for every A G (Aat, A*), and in particular 

u* <u in B. □ 

The following lemma is the key for the proof of the singularity of u* in higher dimensions. 

Lemma 4.2. Let N > 9. Suppose there exist A' > 0, /3 > and a singular radial function w E {B)with 
e LTociB \ {0}) such that 

[ A^w < forO<r<l, 



and 

2/3 



< f (A0)2 for all ^eH^iB), (12) 



1- IfP> A', then X* < A'. 

2. If either /? > A' or if (3 = \' = then the extremal solution u* is necessarily singular. 

Proof: 1) First, note that ^ and e L^ci^ \ {0}) yield to jjt^ e L^{B). By a density argument, 
(fTTj) implies now that u; is a H^{B)—weak sub-solution of {P)x' whenever > 4. If now A' < A*, then 
by Lemma 12.31 w would necessarily be below the minimal solution ux' , which is a contradiction since w is 
singular while ux' is regular. 

2) Suppose first that (3 — X' = and that iV > 9. Since by part 1) we have A* < we get from Lemma 
14.11 and the improved Hardy- Rellich inequality ([8]) that there exists C > so that for all (j) G Hq{B) 



B 



It follows that > and u* must therefore be singular since otherwise, one could use the Implicit 

Function Theorem to continue the minimal branch beyond A* . 
Suppose now that (3 > X' , and let ^ < 7 < 1 in such a way that 

:= (2^)1/3 < 1. (13) 
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Setting w := 1 — a{l ~ w) , we claim that 

u* <w in B. (14) 

Note that by the choice of a we have a'^A' < A*, and therefore to prove (HH) it suffices to show that for 
a^X' < A < A*, we have u\ < w in B. Indeed, fix such A and note that 

» 9 A 9 ctX' a^X' X 
A^iu = aA^w < = -, TTT < 



(l-w)2 (l-iZ;)2 - (l-w)2- 
Assume that ux < w does not hold in i?, and consider 

Ri := sup{0 < i? < 1 I ux{R) > w{R)} > 0. 

Since iu{l) = 1 - a > = wa(1), we then have Ri < 1, ux{Ri) = w{Ri) and (uA)'(i?i) < (iD)'(i?i). 
Introduce, as in the proof of Theorem 13.21 the functions u\ji-^ and wn^^. We have that UA,fli is a classical 
solution of (P)A,a',/3', where 

a' := R;-'iux{Ri) - 1) + 1 , /3' := R^HuxYiRi). 

Since A < A* and then 

2A 2A* 2A' 2/3 

< ^77^ 7^ = ^^ < 



(1 — w)^ a^{l — w)^ 7(1 — w)"^ (l — w)^^ 

by p2)) w/jj is a stable i?2(5)„.^^ga]j sub-solution of {P)x,a',p'- By Lemma [2.31 we deduce that ux > w 
in Bfi^ which is impossible, since id is singular while ux is regular. Note that, arguing as in the proof of 
Theorem 13.21 (a',/?') is an admissible pair. This establishes claim (fT4|) which, combined with the above 
inequality, yields 

2A* 2A* 23 

< ^TT, < 



(1-M*)3 - a3(l-ii;)3 (l-u;)3' 
and therefore 

r (A0)2 _ 2AV 

0effJ(B) ?!)2 

It follows that again fJ.i{u*) > and u* must be singular, since otherwise, one could use the Implicit Function 
Theorem to continue the minimal branch beyond A* . 
Consider for any m > the following function: 

wm:=l- ,^r'" + (15) 
3m — 4 Am — 4 

which satisfies the right boundary conditions: u'm(l) — w'^{l) — 0. We can now prove that the extremal 
solution is singular for > 9. 

Theorem 4.1. Let N > 9. The following upper bounds on X* hold: 

L IfN>31, then Lemma \4. Up holds with w := W2, X' = 27 Xm and (3 = and therefore X*{N) < 27 Xn- 
2. If 17 <N < 30, then LemmaU^holds with w := W3, X' = f3 = ^f, and therefore X*{N) < 



3. If 10 < N < 16, then Lemma \4-2\ holds with w := w^, A^ < /3jv given in Table[^ and therefore 
A* (A) < A'^. 

4-. If N = 9, then Lemma \4.2\ holds with w :— W2.87 ^'9 ■— 366 < /3g := 368.5, and therefore A* (9) < 366. 
The extremal solution is therefore singular for dimension N >9. 
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Table 1: Summary 



9 


W2.S 


366 


366.5 


1 n 
iU 


W3 


45U 


45 ( 


11 


W3 


560 


739 


12 


W3 


680 


1071 


13 


W3 


802 


1495 


14 


W3 


940 


2026 


15 


W3 


1100 


2678 


16 


W3 


1260 


3469 


17 < iV < 30 


W3 


Hn/2 


Hn/2 


iV > 31 


W2 


27Xn 


Hn/2 



Proof. 1) Assume first that N > 31, then 27A < We shaU show that W2 is a singular /f^(_B)— weak 

sub-solution of {P)27\ so that p2)) holds with (3 = Indeed, write 

W2.= l- \x\^ - 2(|a;|3 - \x\^) ^u-cjio, 
where := 2(|a;|^ - |a;p), and note that W2 e H^iB), G L^iB), < W2 < 1 in B, and 

3A ^ 27A 

ri ~ (1 — W2) 

So W2 is /f^(_B)— weak sub-solution of {P)27\- Moreover, by 0o > and ([5]) we get that 



A^W2 = — < 7^— 7T2 in ^ \ {0}- 



for all (/i e H^{B). It follows from Lemma li?^ that u* is singular and that A* < 27A < 
2) Assume 17 < iV < 30 and consider the function 

9 4 43 
:= 1 - -r^ + -r . 
5 5 

We show that is a semi-stable singular iJ^(_B)— weak sub-solution of (P)£jv • Indeed, we clearly have that 

2 

< W3 < I in B, W3 € Hq{B) and G L^iB). To show the stability condition, we consider G Hq{B) 

and write 



-^JV / 7 = 125iJAr / < 125i7Ar SUp 



1 



{1~W3)^ 7b (9ri - 4r3)3 o<r<i (9 - 4ri)3 7s 



a2 



JB ^ Jb 

by virtue of ([8]). An easy computation shows that 

Hn , 2 25i7Ar 9A 12 -1 

A W3 ^ 



2(1-^3)2 2(9rl-4r3)2 5rf 



25iV2(Ar-4)2 8(iV_|)(iV--|) 12 
32(9rl - 4r3)2 5rt 5 
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By using Maple one can verify that this final quantity is nonnegative on (0, 1) whenever 17 < TV < 30, and 
hence is a H^{B)—weak sub-solution of {P)hj±- It follows from Lemma [4.21 that u* is singular and that 



3) Assume 10 < < 16. We shall prove that again w :— satisfies the assumptions of Lemma 14. 21 Indeed, 
using Maple, we show that for each dimension 10 < iV < 16, inequality ^ holds with given by Table [T] 
Then, by using Maple again, we show that for each dimension 10 < A < 16, the following inequality holds 

(Af-2)2(A-4)2 1 (A-l)(A-4)2 1 



16 (|a;|2 - |a;p+i)(|x|2 - 4 |a;|2(|x|2 - ) 

2Pn 



> 



where f3N is again given by Table [T] The above inequality and the Hardy-Rellich inequality ^ guarantee 
that the stability condition (fT2| holds with (3 := (3n- Since f3N > A^, we deduce from Lemma [4.21 that the 
extremal solution is singular for 10 < A^ < 16. 

4) Suppose now A^ = 9 and consider w :— W2.8- Using Maple on can see that 

. 9 366 

^ w < m B 

(1 — ■w)'' 

and 

723 / A^ — 1 \ 

j^^^ < Qir) i^P{r) + -^j for all r G (0, 1), 

where P and Q are given in ifTO]) . Since 723 > 2 x 366, by Lemma [4.21 the extremal solution u* is singular 



in dimension N — 9. □ 



5 Appendix: Improved Hardy-Rellich Inequalities 

We now prove the improved Hardy-Rellich inequalities used in section 4. They rely on the results of 
Ghoussoub-Moradifam in [10] which provide necessary and sufficient conditions for such inequalities to hold. 
At the heart of this characterization is the following notion of a Bessel pair of functions. 

Definition 5.1. Assume that B is a ball of radius R in M^, V,W G C^(0,1), and -pr-^Tyjy^dr = -|-oo. 
Say that the couple {V, W) is a Bessel pair on (0, R) if the ordinary differential equation 

(Bv.w) y"{r) + + ^^)y'{r) + ^y(r) = 

has a positive solution on the interval (0, i?). 

The space of radial functions in (B) will be denoted by (B) . The needed inequalities will follow 
from the following result. 

Theorem 5.1. (Ghoussoub-Moradifam |10| ) Let V and W be positive radial -functions on B\{0}, 
where B is a ball centered at zero with radius R in (N > I) such that ipr^Ty-^dr = +oo and 

r^^^V{r)dr < +oo. The following statements are then equivalent: 

1. (y, W) is a Bessel pair on (0, R). 

2. /^F(|x|)|V0|2dx > ^^W{\x\)(l?dx for all(j)eC^{B). 

3. //linir^o r'^V{r) — for some a < N — 2, then the above are equivalent to 

f V{\x\)iA(f>fdx> f W{\x\)\V(^\^dx + {N-l) f (lM_]iM)|v</,|2dx 
Jb Jb Jb fI FI 

for all4>(^C^^{B). 
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4- If in addition, W{r) — ^^2'''* + — Vrr{r) > on (0, i?), then the above are equivalent to 

v{\x\){Aq)fdx> [ w{\x\)\\/^\^dx + {N -1) [ {YMl^Y:::^)\s/(f,\2dx 



for all 4, e C^{B). 
We shall now deduce the following corollary. 

Corollary 5.1. Let N > 5 and B be the unit ball in . Then the following improved Hardy- Rellich 
inequality holds for all (p G C^{B): 



16 J B {\x\'^ ~\x\-2+^){\x\'^ -\x\-2) 

{n-\){n-aY r (j,^ 

Jb |a;|2(|a;|2-|x|f) 



4 

Proof. Let < a < 1 and define y{r) := — a. Since 

y" + _{N- 2)2 1 



y 

the couple ( 1, ^^^^^^ ^v — - ) is a Bessel pair on (0, 1). By Theoreni l5.1f 4) the following inequality then 

holds: 

4 Ji3 - a\x\^+^ Jb \A 

for aU (f> e C(f (B). Set F(r) := and note that 

V r 2 1 _ ar^~^ ~ f 

The function y(r) = r^-T+'^ — 1 is decreasing and is then a positive super-solution on (0, 1) for the ODE 

// ,N^1 Vr, , Wi{r) 

where 

W^ M = — 

Hence, by Theorem 15.1^ 2) we deduce 



If +1 - ^ 2 



B {\x\^ ~ a\x\^+^){\x\^ - |x|f ) 



for all (f) G C^{B). Similarly, for V{r) = we have that 

12 — o 



B - ^ 2 '7s |x|2(|a;|2 - |x|f ) 

for all 4) G C^{B). Combining the above two inequalities with p?|) and letting a 1 we get inequality 
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Corollary 5.2. Let N = 9 and B be the unit ball in R^. Define (p{r) :— r + r — 1.9 and ip{r) := 

r^"2"+^ + 20r^^'^^ + lOr^^ + lOr + 7r^ — 48. Then the following improved Hardy- Rellich inequality holds for 
all(j)(zC^{B): 

{Acj^f > Qi\x\) (^Pi\x\) + ^) (18) 
where 

^"ir) + E_l^'{r) r{r) + ^j^' {r) 
P{r) and Q{r) := — ^ . 

Proof. By definition (l,P(r)) is a Bessel pair on (0,1). One can easily see that P{r) > Hence, by 
Theorem 15. If 4) the following inequality holds: 

(A0)2dx> / P(|a;|)|V0|2 + (7V-l) / (19) 

for all (f) G (B) . Using Maple it is easy to see that 

^>-- in (0,1), 
P r 

and therefore ip{r) is a positive super-solution for the ODE 

A^-1 P,(r) F(r)Q(r) 
r " (^) "('') 

on (0, 1). Hence, by Theorem O; 2) we have for all (j> £ C^{B) 

P{\x\)m'> [ P{\x\)Q{\x\)cp', 

JB JB 

and similarly 



B - Jb 



since ip{r) is a positive solution for the ODE 

y" + ^^y'{r)+Q{r)y^O. 
r 

Combining the above two inequalities with (fT9|l we get (fT8|l . □ 
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